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We show that three distant labs A, B and C, having no prior entanglement can establish a shared GHZ state, 
by mediating two particles which remain separable from each other and from all the other parties throughout 
the process. The success probability is 1. We prove this in a general framework for systematic distribution of 
entangled states between 2 and more parties with separable states in d dimensions. The proposed method may 
facilitate the construction of multi-node quantum networks and many other processes which use multi-partite 
entangled states. 


Introduction Decades of study of entanglement, from 
its inception in relation with conceptual framework of the 
quantum theory [1-3], to its recent upsurge as a tool in 
quantum computation and quantum information science [4] 
and the inevitable need for its quantification has shown that 
its various unexplored aspects are still intriguing and can 
still surprise us. In the history of study of entanglement, 
the turning point was the discovery that entanglement can 
be used to teleport quantum states [5], to share secret keys 
and to do dense coding, all pointing to the fact that this 
concept is a useful resource, as important as many other 
physical quantities. Since then it has been shown that entan¬ 
glement can be manipulated [6], measured [7], distributed[8] 
and distilled [9]. The distinctive feature of entanglement 
is that it cannot be created by local action and classical 
communication (LOCC) [10]. Any attempt for distributing 
entanglement between two or more distant particles, should 
involve either a direct interaction between the particles when 
they are close [11] or should be mediated between them 
by sending a mediating particle through a quantum channel 
[12]. This latter method is of course very much vulnerable 
to noise, since the state of mediating particle is very frag¬ 
ile. This problem becomes very severe if we note that the 
mediating particle may be entangled with the distant particles. 


It came as a big surprise when it was shown [13] that a 
mediating particle (C) can be entirely separable from the 
distant particles (A and B) which are to be entangled. It 
was then shown that such a protocol is also possible with 
Gaussian states [14, 15] where the actual experiments with 
these continuous variable states were reported in [16, 17] and 
for qubit states of light in [18]. The question of what type of 
initial qubit states can lead to this phenomena was explored 
in [19] and the cost for this entanglement increase was related 
to the discord [20, 21] in the initial state in [22, 23]. 

However important questions have remained unanswered 
until now: Can we distribute multi-partite entanglement 
between N parties by using separable ancillas? If yes, how? 
How do we distribute d— level maximally entangled states 
by using separable ancillas? Can we provide a systematic 
method for these protocols? These questions are not only 
of conceptual relevance, specially in view of the fact that 
GHZ states [24] have a different kind of non-locality than 


Bell states, but they have practical relevance in view of 
their role in many quantum protocols like quantum secret 
sharing [25, 26] and other quantum protocols using quantum 
networks [27]. Moreover in view of the studies in [22] and 
[23] on the cost of entanglement sharing and its relation with 
quantum correlations, a systematic generalization of this 
phenomena to d levels and N parties, will certainly lead to 
deeper insight on these inter-relations and the bounds that 
they suggest. In this letter we provide positive answers to all 
these questions. We first discuss d— level Bell state distribu¬ 
tion and then go on to GHZ distribution between many parties. 

Notations and conventions: 

• In all the multi-partite states, we write the order of 
parties only in the left hand side of the states which 
is meant to be the same as in the right hand side. For 
example we write the same state as \^)abc = |011) or 
I ^)bca = 1110), since by necessity we have to change 
the order of terms in the states, the reader is asked to 
always check the order of party labels on the left hand 
side. A CNOT operation on C, controlled by A is 
denoted by C AC and acts as C AC \i, j) AC = |i,i + j). 
The combined states before Alice CNOT operation(s), 
after her CNOT operation(s) and after Bob CNOT 
operations are denoted respectively by p(°\ p ^ and 


Distribution of Bell states in arbitrary dimensions Con¬ 
sider two parties, Alice(A) and Bob(B), who possess pure 
states of the form 

i mu = 4jX> Sife ii>, i <K-k)) B * ''i> Sjk \j)- 

Vd J=0 v d J=0 

(1) 

where K is an integer, cc = is a phase, 

k = 0,1,•••AT — 1, and the non-negative integers Sj 
are all different from each other and satisfy the following 
conditions: 

i) 0 < Si < y 
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ii) for any four indices i, j, m and n, which are not all equal 

Si+Sj = s m +s n , only if (i,j) = ( m,n ) or (i,j) = ( n,m ). 

( 2 ) 

These conditions will play an important role in the process as 
we will see. In particular we note the consequences 


CNOT operation). If we disentangle C in each of these two 
states, it means that C has been a separable ancilla in its 
travel between Alice and Bob during the whole process. 

Consider first the states |T>^^(/c)) in (3) and determine the 
un-normalized mixture, 


1 K ~ x 

K ^ h3S 


k =0 
K-l 


1 

— LJ ( S i + Sj- 8 m-s n )k _ x. X. , X X. _ X. .£ X. 

— u i,m u j,n \ u i,n u j,m u ij u mn u jm 


k=0 


Alice has also a d-level ancilla C at state |0)c- She now 
enacts the operation C AC , turning the fully separable state 
I^))JO)cI<K-*0)b i nt0 the following partially entangled 
state, where C is entangled to A but not to B: 


l$ (1) (fc)),OB = ^ E (3) 

j=o,i=o 

and sends the ancilla C to B where he performs C~* , turn¬ 
ing the state to 

i^ (2) (fc)>c^ = \ E “ {sj ~ si)k \j-ij,i) 

j=o,i=o 

= -4 EHbU*0> W 

m=0 

where 

\x m (k)) AB = -m), (5) 

d 3=0 

are maximally entangled states. The important point, as we 
will see, is that |x 0 ) and only this state, is independent of k. 

If now the ancilla qudit C is measured in the computational 
basis, the state of AB is projected onto the maximally 
entangled state \x 0 )ab with a given probability. In this way 
A and B, have used the ancilla C to entangle their initially 
separable state \(/>(k))k))s into a maximally entangled 
one. The problem however is that the ancilla has been 
entangled to the states during the process. We now show that 
this entanglement can entirely be removed and the qudit C 
can be completely separated from the state of the two parties 
AB in all stages of the process. 

To do this we use the fact that a mixture of entangled 
states can be separable and use two consecutive processes, 
namely mixing and symmetrization. The two main stages of 
the process correspond to the two states \^ <yl \k)) ACB ( after 
Alice CNOT operation) and \^ 2 \k)) CAB (after Bob inverse 


K -1 

pTd = E \* {1) m(* {1) m (6) 

k =0 

which by using the conditions on becomes 


K 

d? 


p'^ B = ^\T.^3i+d\GHZ){GHZ\) 

/ 


where 7= | ijk)(ijk\ and 

\ghz) acb = 2= E \ jjj)- 


(7) 


( 8 ) 


In view of the GHZ term the ancilla C has not been made 
separable from AB. However we now note that by the very 
operation of Alice on C, which involves only A and C, this 
state has the structure of p AC , that is B is certainly separa¬ 
ble from AC. Therefore if we make it symmetric with respect 
to the interchange B ^ C, then it will have also the structure 
AB\C and the ancilla C will be separable from the original 
particles A and B. The GHZ state has obviously this symme¬ 
try, so we add the projectors i\jij to this state and obtain the 
separable state 



which after proper normalization gives us the separable den¬ 
sity matrix after operation of Alice and before Bob operation 

A'L = + 

( 10 ) 

From this state and by noting that it can be written with both 
forms of indices ACB or ABC (due to its symmetry) and 
by inverting the CNOT operation of Alice, we will obtain the 
fully separable original state of Alice, Bob and the ancilla: 

= K( J_Yj E W*»W*)l <8 w-fc»w-t)l 810){0| 
+ ( 11 ) 

From (10) we now obtain the final state after Bob’s inverse 
CNOT operation, 

= 25 ^Tlxo><xol®|o><o| 
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+ d(2d- 1) + ( 12 ) 

This shows that when the ancilla is measured in the compu¬ 
tational basis, the state of AB is projected onto a Bell state 
with probability 7 ^—^. This probability is independent of the 
number of states K. This number is determined once we have 
chosen the parameters Sj satisfying the two conditions. We 
illustrate this with a few examples: 

Example 1- qubits : The simplest assignment is to take 
80 = 0 and si = 1. With K = 3, we obtain uj = e~^, 
where Alice and Bob each use three states. Another solu¬ 
tion is given if we take K = 4 or cj = e 22 ?, where Alice 
and Bob each use 4 states. This leads to the state given in [13]. 

Example 2- qutrits : The simplest assignment is 

so = 0,8i = 1,82 = 3. The smallest number of states 
is K = 7 with uj = e "T 1 . 

Example 3- four-level states : The simplest assignment is 
80 = 0,8i = 1,82 = 3, 83 = 7. The smallest number of 
states is AT = 15. 



FIG. 1: (Color online) GHZ distribution with separable states, (a) 
Alice operates on two ancillas with her CNOT operations, controlled 
by her qubit A and sends the two ancillas to Bob and Charlie, (b) they 
operate on their respective ancillas by their CNOT’s (inverse CNOT’s 
in case of qudits). (c) At the end a GHZ state is shared between A, 
B and C with a probability Psuccess — \ • Throughout the process 
the ancillas remain separable from each other and from ABC. The 
process can be generalized for N parties and for qudits. 


Example 4- d level states : The simplest assignment is given 
by Si = 2 l — 1 and K = 2 d — 1 states. 

Distribution of GHZ states We now use the same strategy 
to distribute GHZ state between three parties, figure (1). The 
idea is again mixing and symmetrization of product states. To 
simplify the presentation we explain in detail the idea for 2- 
level GHZ states between 3 parties. This simple setting allows 
us to understand where the conditions on the state parameters 
come from and how one should proceed for other cases. In 
view of the previous section, generalization to d level states 
and more parties is straightforward, although its presentation 
will clutter the text with many indices. 

Here Alice, Bob and Charlie possess pure states of the form 
^(|0) + OL k |1», ^(|0) + /3fc|l))j and ^=(|0) + 7fc|l)) re¬ 
spectively, where Pklk — 1- Alice has two ancillas Di and 
D 2 which are both at the states |0). Alice now performs the 
operations C AL>1 and C AE>2 and turn the product state into 

I $ (1) (fc)U lD2 ,BC = 

-F(|000) + a fe |lll»(|0) + /3 fe |l»(|0) + 7fc|l»-U3) 

and then sends the ancilla qubits D\ and D 2 to Bob and Char¬ 
lie who perform their respective CNOT operations (inverse 
CNOT for d— level states) and turn the state into 

D lD2 ,ABC = ^(|00)|Xoo( fc )) + l 01 )IXoi( fc )) 

+ |10)|x 1 o(fc)) + |ll)lx 11 (fc))) (14) 

Ixoo) = A(|000> + |111», 


lx°i) = 2^(7|001) +7- 1 |110)), 
lx 10 ) = A^(/3|010) +/3- 1 |101)), 

Ixii) = A( a |ioo) + a -i|oii)). 

Note that |xoo) is the standard GHZ 3 state and it is important 
that it is independent of k. 

Now Bob and Charlie can measure the ancillas D\ and D 2 
in the computational basis which projects their joint state with 
A into the GHZ state |x 00 ) with a given probability (to be 
determined later). To remove the entanglement of the two an¬ 
cillas Di and D 2 from the rest of states we again use the pro¬ 
cedure of mixing and symmetrization. We first mix the final 
state (14) and form p'^ D ^ ABC := |$ (2) (*0>( < l ,(2) (fc)l- To 

make this state separable we require that 

Z \Xij( k ))(Xi’j'(k)\ = 0 (i,j) ^ ( i',j '), (15) 

k 

which leads to the independent equations: 

b b b b b 1 ^ b ' ^ 

r\j r\j r\j r\j rv 

(16) 

Under these conditions, the mixed state will be 

pff lABC = ||00)(00| ® \GHZ 3 )(GHZ 3 \. (17) 

We will later come back to this state. The next step is to 
mix the states (13) after Alice’s and before Bob’s opera¬ 
tion. So we form the un-normalized mixture p 'W := 

r AD 1 D 2 ,BC 


where 
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J 2 k |4>( 1 )(A;))(<I>( 1 )(A;)| and symmetrize it with respect to the 
interchange (DiD 2 BC). Since by the very operation 
of Alice on the ancillas alone, this state has the structure 
Pad 1 d 2 \bc ’ if we can ma k e it symmetric, it will have also the 
structure of p ABC \ DlD2 , hence the ancilla will be separable 
from ABC. Using an abbreviated notation as V = D\D 2 and 
B = BC and also using the binary notations for two-qubit 
states (e.g. |3) = 111)) we can write 

\® W (k)) Ave = ^(|000)+ 7 |00T)+/3|002)+/3 7 |003) 

+a|130) + cry |13l) + a/?|132) + 1133)). (18) 

In view of this new form of the state |3>( 1 )(fc)), we note that 
the non-symmetric terms in p'^ are of two kinds: 1) those 
which are projectors like 7r 0 Q2 and 2) those which are cross 
terms like 1000) (0021. We demand that all these cross terms 
vanish. A glimpse at (18) shows that this requires that 

£“/= = £/3i = £ 7 ;? = o, (!9) 

k k k 

where we have also used condition (16). We are then left with 

3 2 

2| GHZ 5 )(GHZ 5 I + £ TTooi + £ TT^) 

i=1 i =0 

which can be symmetrized by adding suitable projectors. Af¬ 
ter normalization the state will be 

p ( 1L = \\ghz 5 )(ghz 5 \ 

3 2 

£(*«H + ^oio) + £(*13? + n ia) • 

,i=l i =0 

Note from this state that the two particles D\ and D 2 , are not 
only separable from A, B and C, but are also separable from 
each other. 

We should now consider what effect this addition of 
projectors has on the initial and p( 2 \ In view of the 
fact that the effect of CNOT’s by Alice and Bob turn any 
projector into another projector, we see that by this addition 
the initial state of Alice, Bob and Charlie and also 
their final states p^ remain separable and the probability of 
distributing a GH Z 3 states is given by P SU ccess = \ • 

It remains to solve the constraints (16) and (19). To do this 
we use the relation OLkfiklk — 1 to express everything in terms 
of otk and (3k and make the ansatz that (3k = a\ to transform 
the constraints (16) and (19) into 

£<4 = o, r = 1 , 2, • • • 6. (20) 

k 



= * 

'ABB O 


This is simply solved by taking ak = e 7 

27 xik Airik _ 

ak = e 7 , p k = e 7 , lk = eJW 


which then gives 

87 rik 


( 21 ) 


and k ranges from 0 to 6 , that is, Alice, Bob and Charlie 
each should use 7 different states of the form given in the 
beginning of this section. 

The same analysis can be carried out for N parties. In fact 
for N parties, each party Aj, j = 1 , 2, • • • N, should have 
2 ^ — 1 states of the form 

\Mk))A, = ^(\ 0 )+^ k \l)), k = 0 ,l, 2 ,--- 2 N - 2 , 

( 22 ) 

2 7T 7 

where uj = e 2N - 1 . It is easy to show that running this pro¬ 
tocol in the same way as before leads to sharing a GHZ state 
between the parties with probability P SU ccess = 2 n_ 1 • 

In summary, we have introduced a constructive method for 
using separable states to distribute d-level Bell states between 
2 parties and GHZ states between 3 and more parties. The 
method is straightforward to generalize to GHZn sharing 
for qudits. Several questions remain for future investigations: 
First one can consider other classes of multi-partite states, 
like W states, 1-uniform states [28] and cluster states [29] 
which are important for one-way quantum computation. 
Second one may see how GHZ state sharing with separable 
states can be implemented with Gaussian variables along the 
the theoretical lines set in [14, 15] and their experimental 
realization in [16-18]. Finally, investigation of the cost of 
multi-partite entanglement sharing [22], [23] and its relation 
with multi-partite quantum correlation present in the initial 
states, will surf again, a study which may lead to a host of 
connections and bounds between two not-well understood 
concepts. 

We thank Shakib Vedaei for interesting discussions. 
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